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Indian monsoon rainfall is higher in epochs of higher solar 
activity: a wavelet cross-spectral analysis 
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Using the Morlet continuous wavelet transform on 
data over the period 1871-1990, it is found that the global 
wavelet cross spectra between two solar activity indices 
and seven major Indian monsoon rainfall time series show 
significant power around the period of the 11 year so- 
lar cycle, passing the test of significance proposed by 
Torrence and Compo (1998) at levels exceeding 95% in 
10 out of the 14 cases studied. Furthermore two dis- 
tinct epochs are found in the computed cross-spectrum, 
the transition between them occurring around the decade 
1915-1925, beyond which solar activity parameters show 
a significant increase. By comparison between selected 
periods of three cycles in solar activity in each of the two 
epochs, it is shown that the average rainfall is higher in 
all seven rainfall indices during periods of greater solar 
activity, at 2;-test confidence levels of 95% or greater in 
three of them. 



1. Introduction 

There are numerous studies of the possible influ- 
ence of solar activity on terrestrial climate. Although 
solar-terrestrial connections have been somewhat scepti- 
cally received by much of the meteorological community, 
the growing availability of paleo-climate indicators us- 
ing novel measurement techniques, and the discovery of 
the prominent solar activity period of 11 years in var- 
ied climate records using new mathematical tools, are 
producing increasingly strong evidence for possible asso- 
ciation between solar processes and terrestrial weather 
and climate indices (Beer et o/.(1990), Friis-Christensen 
and Lassen (1991), Lassen and Friis-Christensen (1995), 
Haigh (1996, 1999, 2001), Labitzke and Loon (1997), 
Mehta and Lau (1997)). Doubts about effects on shorter 
time scales (of the order of lOOy) seem however to re- 
main. Earlier studies of possible connections between 
Indian monsoon rainfall and solar activity (Jagannathan 
and Bhalme (1973) and Jagannathan and Parthasarathy 
(1973)), using correlation and power spectral analysis 
of the rainfall distribution in 48 meteorological stations 
spread all over India, reported presence of the 11-year 
sunspot cycle at signiflcance levels of 95% or higher in 5 
of them. 

The possibility of shedding light on the variability of 
the Indian monsoon using wavelet techniques (Kailas and 
Narasimha (2000), Narasimha and Kailas (2001)), and 
the idea that the tropics can amplify a small radiant flux 
signal to a relatively large and dynamic climate change 
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elsewhere in the world as well (Haigh (2001), Visser et 
al. (2003)), serve to provide further motivation for the 
present work. In a recent study by Neff et al. (2001), 
strong coherence between the solar variability and the 
monsoon in Oman in the period between 9 and 6 kyr 
B.P. has been reported. 

There has been a recent surge in interest in incorpo- 
rating wavelet techniques for climate signal detection and 
analysis, as can be seen in the works of Lau and Wcng 
(1995) and Torrencc and Compo (1998). In particular, 
Narasirnha and Kailas (2001), analysing Indian monsoon 
rainfall, identified six 'qucisi cycles', with average periods 
of 2.92, 5.8, 11.4, 19.9, 34.2 and 79.8 y. These num- 
bers suggested the presence of several super- and sub- 
harmonics of the basic sun-spot period of about 11 y. 
Iyengar and Raghukantha (2004) have identified intrinsic 
modes with periods virtually identical to those observed 
by Narasirnha and Kailas (2001). 

In contirmation of this on-going research, we now 
present an analysis of the association between indicators 
of solar activity and Indian monsoon rainfall, using the 
continuous wavelet transform method with the Morlet 
wavelet function. 

2. The data analysed 

Seven annual area-weighted rainfall time series for the 
period 1871-1990 have been considered for the analysis, 
namely, all India summer monsoon rainfall (AISMR), 

north east India (NEI), north west India (NWI), cen- 
tral north cast India (CNEI), west central India (WCI), 
peninsular India (PENSI) and the homogeneous Indian 
monsoon (HIM)(Parthasarathyet al (1995)). The choice 
of these regions is based on several well defined crite- 
ria including primarily a considerable degree of spatial 
coherence. The HIM region covers the central and north- 
western parts of India amounting to 55% of the total land 
area of the country, and may be seen as the most char- 
acteristic index of the component of Indian rainfall dom- 
inated by the dynamics of the south west monsoon. As 
it is well known that the NEI region shows in many ways 
an anti-correlation with the HIM region (Parthasarathy 
et al. 1993), we present more detailed results specifically 
for these two regions. 

The range of scales over which these rainfall data can 
provide useful information on temporal variability is lim- 
ited at one end by resolution, since not more than 12 
points per year are available, and at the other end by 
the limited length of data stretch of 120 years. For the 
present study wc have found annual rainfall to be the 
most appropriate rainfall index to use. The solar in- 
dices under study arc sunspot number index and group 
sunspot number. Sunspot areas have also been stud- 
ied, but have been found to provide no new informal 
tion. The sunspot index data have been obtained from 
Rai Choudhuri (1999) and Fligge et al. (1999), and the 
data for group sunspot number from the NOAA ftp site 
ftp://ftp.ngdc. noaa.gov/STP/SOLAR.DA TA. 

We must note one important difference between the 
rainfall and sunspot data. The former are cumula- 
tive, the monthly data being a sum of the daily data. 
That is not the case for sunspot data, since individual 
sunspots live typically for several days. All monthly 
data of sunspots are usually compiled by taking av- 
erages of daily data over a month. The monthly av- 
erage sunspot number plotted against time does not 
appear very smooth (see, for example, the website 
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http: / /science. msfc.nasa.gov /ssl/pad / solar / image /zurich.gif ) . 
As our aim here is to study possible correlations of the 
rainfall data with solar processes with time scales of order 
years to decades, we filter out fluctuations in the sunspot 
data at small time scales by using yearly sunspot and 
group sunspot numbers. 



3. Wavelet Cross Power Spectrum 

We use the Morlet wavelet function 

V)(r)) = 7r-l/*e*'^'"'e-'''/^ (1) 

where too is a nondimensional frequency, taken equal to 
6 in order to satisfy the wavelet admissibility condition, 
and r? is a nondimensional time parameter. For a discrete 
sequence i?„, n = 0, .., {N — 1), the continuous wavelet 
transform W^{s) is defined as a convolution of Rn with 
a scaled and translated version of the wavelet function 
ip{ri), as given by the expression 

w^.%)=i:^„'V'n^^^^], (2) 



where * denotes the complex conjugate, and 5t is the 
(sampling) time interval between two consecutive points 
in the time series. The wavelet function at each scale s 
is normalised to have unit energy, so that 

tpisuJk) = (^) V'(s'^fc), (3) 

I i,{co') \^ dJ = 1, (4) 

where 'ip is the Fourier transform of The wavelet power 
spectrum of Rn is given by the convolution Wjf (s) * 

Wn^"* and the wavelet power is given by the magni- 
tude 1 W"{s) * [W;f (s)]* I- For the wavelet cross power 
spectral analysis wc utilize the casy-to-usc toolkit, includ- 
ing statistical significance testing, as outlined by Torrence 
and Compo (1998). The cross wavelet spectrum W^^ {s) 
between two time series Rnit) and Snit), with the re- 
spective wavelet transforms Wn{s) and Wn{s), may be 
defined as 

W^%s) = W^{s)[wUs)Y . (5) 

The cross wavelet power is | Wn^{s) \. Torrence and 
Compo (1998) derive the confidence levels for the cross 
wavelet power from the square root of the product of two 
distributions. In the test, a peak in the wavelet power 
spectrum is considered to be a true feature, with a certain 
percentage confidence, if the peak is significantly above 
a background or reference spectrum given by 

^ l-hQ2-2acos(27rfc/iV)' ^ ' 

Here k is the frequency index, and a = {a\ + y^2)/2 
where qi and 02 are the lag-1 and lag-2 autocorrela- 
tion coefficients of the process under consideration. For 
a white noise background spectrum a = 0. If the two 
time series have background spectra given respectively 
by and , then the cross wavelet power distribution 
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will be given by 




{s)\ 



(7) 



where aR and as are the standard deviations of R and S 
respectively, v is the number of degrees of freedom with 
which is distributed, p denotes the level of confidence 

and Z„{p) denotes the value of the distribution with 
v degrees of freedom at the confidence level p. For the 
complex Morlet wavelet u = 2 . 

4. Results 

A plot of the annual time scries of the four rainfall and 
two solar indices under consideration is shown in figure 
1. The rainfall time scries appear irregular and random, 
while the solar indicators have a clearly cyclic character. 

The results of the present analysis are presented in 
the form of colour-coded contour maps (see figures 2 and 
3) of wavelet cross power spectra as functions of time 
and Fourier period (henceforth refered to as period) re- 
spectively for HIM and NEI rainfall. Outlined on these 
graphs arc thick contours enclosing regions where wavelet 
cross power exceeds 95% confidence levels, with respect 
to the reference spectra mentioned above. The cones of 
influence within which edge effects become important are 
also indicated by dashed lines in the figures. 

Figures 4 and 5 show several global wavelet power 
spectra as functions of period for HIM and NEI rain- 
fall. For HIM rainfall (figure 2), a wavelet cross power of 
noticeably high magnitude is observed at the middle of 
the period range 8-16 years. An integration over time of 
this wavelet cross power, which gives the global wavelet 
cross power spectrum, is shown in figure 4. It can be 
seen clearly from this figure that the cross wavelet power 
crosses the 95% confidence line at a period of around 11 
years. Circles and squares respectively denote the con- 
tributions of the individual global wavelet power of the 
rainfall and the sunspot time series to the global spec- 
trum. It is seen that the contribution of the 8-16 year 
period band is 16.8% for the total rainfall and 82.8% for 
the sunspot index. 

Figures 3 and 5 show plots similar respectively to fig- 
ures 2 and 4 but for NEI rainfall and group sunspot 
number. Significant power is observed in the 8-16 year 
period band in this case also. The global cross wavelet 
power spectrum also crosses the 95% confidence line in 
the period range 8-16 years (figure 5); the correspond- 
ing cross power (see figure 2) is however lower than that 
for HIM. The highest cross power is observed during pe- 
riods 1890-1905 and 1915-1965 in the former case, and 
during 1880-1920 and 1935-1980 in the latter case. In 
both cases, therefore, there is strong indication of two 
relatively distinct epochs, the transition from one to the 
other occurring around 1915-1925. As may be seen from 
Figure 1, solar activity generally shows an increase be- 
yond this decade, suggesting that higher solar activity is 
associated with the higher cross spectrum. 

Similar analyses have been carried out for the other 5 
rainfall indices as well, and the results are summarised in 
Table 1. It wiU be seen that, except for NEI and PENSI 
rainfall, all the other cases show up regions of cross power 
at confidence levels exceeding 95% at the 8-16 year period 
band. 

From figure 1, two epochs of low and high group 
sunspot number can be identified, this transition also oc- 
curring around 1915 to 1925. We select here one time 
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interval in each of the two epochs, namely 1878-1913 and 
1933-1964 respectively, during each of which three com- 
plete solar cycles are present (considering the period from 
one minimum to the next as one complete cycle). Table 
1 lists the annual rainfall means fii and fj,2 respectively 
for the two periods 1878-1913 and 1933-1964; both the 
epochs and the corresponding means are shown in figure 
1. The null hypothesis that the difference in mean annual 
rainfall between these two periods is zero is rejected at 
the maximum confidence levels listed in Table 1 using a 
onc-tailcd z-test (Crow et al. 1960). Thus the mean an- 
nual rainfall during 1933-64 (higher solar activity) is ev- 
erywhere higher than that during the period (1878-1913) 
of lower solar activity. However the confidence levels are 
95% or higher in 3 cases out 7, including AISMR and 
HIM, and reaching 99% in WCI. At the other extreme, 
it is a low 75% in NEI and NWI. 

The ratio of wavelet power present in the period 8-16 
years to the total rainfall wavelet power is also presented 
in table 1. So arc the ratios of cross wavelet power be- 
tween the rainfall and group sunspot number in the 8-16 
year band to the total cross wavelet power. In the case 
of HIM and WCI, which show the solar-monsoon link at 
the highest confidence levels, the cross spectrum is also 
high (about 56%), and the contribution from the 8-16 y 
band is nearly 17%. It will be seen from Table 1 that, in 
general, significant effects on rainfall go with significant 
levels of cross power. The exception is the arid region 
NWI, where the lower significance levels for the differ- 
ences in rainfall are due to the large relative standard 
deviation. 

Incidentally the present results demonstrate the ad- 
vantages of the wavelet approach, as compared with 
classical correlation/power spectral density methods: (i) 
wavelets permit identification of epochs during which cor- 
relations at different significance levels may have pre- 
vailed; (ii) wavelet methods allow us to take account of 
slight variations iu the effective period or scale ('mean- 
dering') of the effect of a given forcing (as may be seen 
from the regions of high cross spectra in figures 2 and 3), 
such meandering being presumably the result of the non- 
linear interactions between different modes of the system. 



5. Conclusions 

The present study, involving two solar index time 
series and seven Indian rainfall time series and using 
wavelet cross spectral density analysis as outlined by Tor- 
rence and Compo (1998), reveals considerable power in 
the global cross power spectra around the 11 year solar 
cycle period for all the indices considered. In particular 
the global cross power spectra for AISMR, WCI and HIM 
rainfall with the group sunspot number reveal a signifi- 
cant peak at the 11 year period at confidence levels of 95% 
or higher. Greater solar activity seems to be associated 
in all cases with greater rainfall, although at significance 
levels that are distinctly high in 3 and lower in 4 out of 7 
cases studied. This regional variation is not inconsistent 
with the simulations of Haigh et al. (2004), which sug- 
gest that a major effect of higher solar activity may be 
a displacement in the Hadley cell. Such a displacement, 
depending on its magnitude, can have different effects on 
rainfall in different regions. 
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Table 1. Confidence levels and % wavelet power 



Region 


mm 


M2 
mm 


% con- 
fidence 


% power 
rainfall 


% cross 
power ° 


% conf., 
cr-power ^ 


AISMR 


853.5 


883.0 


95 


10.2 


49.9 


98 


HIM 


858.6 


916.5 


96.8 


16.8 


56.7 


99 


WCI 


1067.2 


1145.7 


99 


16.9 


56.4 


98 


PENSI 


1140.5 


1183.2 


85 


10.4 


47.4 


92.5 


CNEI 


1204.2 


1235.4 


80 


7.7 


45.7 


95 


NWI 


542.0 


565.1 


75 


15.3 


55.9 


99 


NEI 


2071.8 


2100.2 


75 


8.6 


45.8 


98 



^ Mean rainfall over three cycles of low solar activity, 1878- 
1913. 

^ Mean rainfall over three cycles of high solaj: activity, 1933- 
1964. 

Confidence level at which m — /i2 differs from zero by the 

z-test. 

Contribution to total rainfall from 8-16 y band. 

" Contribution to cross-spectrum power from 8-16 y band. 

' Approximate maximum confidence level at which cross 
spectrum exceeds values for reference spectrum over a contin- 
uous period of at least 10 y. 
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Figure 1. Raw time series of rainfall and solar indices, 
indicating epochs of low and high solar activity selected 
for analysis, and the means over the epochs for each pa- 
rameter. 
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Figure 2. Wavelet cross power spectrum between HIM 
rainfall and group sunspot number, contour at 99% con- 
fidence level 
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Figure 3. Wavelet cross power spectrum between NEI 
rainfall and group sunspot number, contour at 90% con- 
fidence level 
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Figure 4. Global wavelet cross power spectrum between 
HIM rainfall and group sunspot number 
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Figure 5. Global wavelet cross power spectrum between 
NEI rainfall and group sunspot number 



